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Abstract 

In the framework of a specific scheme of the QCD sum rules for S-wave levels 
of the heavy quarkonium, one derives expressions, relating the leptonic constants, 
£f~) | the energetic density of quarkonium states and universal characteristics in the 

heavy quarkonium physics, such as the difference between the masses of heavy 
quark Q and meson (Qq) and the number of heavy quarkonium levels below the 
^ threshold of (QQ) — ► {Qq) + {Qq) decay. 

in 

^ ■ Introduction 

a: 

Ph! Powerful tools in studies of heavy quarkonia, the bound states of two heavy quarks, are 

r-| ! phenomenological potential models ^ ||, [| and QCD sum rules An applicability 

of the approaches to the systems of two heavy quarks is caused by 

1) a low value of the ratio Kqcd/^q <C 1, where tuq is the heavy quark mass and 
A.qcd is the quark confinement scale, determining the inverse distance between the 
quarks in the bound states, and 

2) a nonrelativistic motion of the heavy quarks inside the quarkonium, v — ► 0. 
In the QCD sum rules, the low value of ratio Aqcd/t^q determines not large contri- 
bution of higher orders of the QCD perturbation theory over the quark-gluon coupling 
as ~ 1/ \n(mQ/ Aqcd) <C 1 in the expansion of Wilson's coefficients, and it makes a 
suppression of nonperturbative quark-gluon condensate contribution, having the power 
form, as 0(AQ CD /mQ) for the gluonic condensate (as G*) contribution into the sum 
rules for vector currents, for example, at some low and moderate values of a number 
for the moments of the spectral density @, f| . 

As has been shown in ref . , the nonperturbative effects in the heavy quarkonium 
spectroscopy can not be described by a potential model, since the correct description 
must take into account the retardation of interaction with condensates. Nevertheless, 
the potential models are suitable for the phenomenological studies, because they can 
reproduce the model formulae or numbers for the quantities, used as input values 
(the level masses, for example). Hence, the potential models can be considered as 
phenomenological meaningful fittings of some experimental values, but they can not 
restore a true potential, that does not exist due to the nonperturbative effects. 
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In the potential models, from data on the spectroscopy of the (cc) charmonium and 
the (66) bottomonium one finds that the nonrelativistic quark motion allows one to get 
the phenomenological potential in the range of average distances between the heavy 
quarks inside the quarkonia 

0.1 fm < r < 1 fm . (1) 

Being the potential of static sources for the gluon field, this potential must not 
depend on the flavours of sources. This flavour-independence is empirically confirmed 
for the QCD motivated potentials [|]]. Such potentials, possessing different asymptotic 
properties in the regions of r — > and r —>■ oo, coincide with each other in the 
region ([!]), where they allow approximations, having a simple scaling behaviour. These 
approximations are the logarithmic || and power || laws 



By the virial theorem 



one finds 



V L (r) = c L + d L \n(A L r) , 
V M (r) = c M + d M (A M r) k . 

/T\ l l dV \ 

(Tl) = di/2 = const. , 
(Tm) = ^(-c M + E), 



(2) 
(3) 

(4) 

(5) 
(6) 



where E is the binding energy of quarks in the quarkonium. Phenomenologically, 
one has k 1, \E\ \cm\, so that in the region of average distances between the 
heavy quarks in the heavy quarkonium ([!]), the kinetic energy of quarks practically is 
a constant value, independent of the quark flavours, 

(Tm) ^ const. (7) 
Then from the Feynman-Hellmann theorem 

f — ®., (8) 

where \i is the reduced mass of heavy quark system (QQ') , one can get that the level 
difference in the system does not depend on the reduced mass of quarks, i.e. on the 
quark flavours, 

E(n, n) — E(n, /i) = E(n, fi') — E(n, u') . (9) 

Condition (|9|) means that the energetic density of heavy quarkonium levels does not 
depend on the quark flavours 

£ = (10) 
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where 4>{n) does not depend on /i. The Bohr-Sommerfeld procedure for the quantiza- 
tion of the nonrelativistic systems with the logarithmic and power potentials gives 



dMh IT 1 dM, 



n 



71=1), (11) 



dn n n dn 

M 

(12) 



dMl = 1 /n 2 /i \ 
dn n V jJt, J 



Since k -C 1, the 5-state density in the Martin potential only logarithmically depends 
on the reduced mass \i. With the same accuracy, one can state that eq.flTTD is approx- 
imately valid in the system with the Martin potential. 

The described properties of heavy quark potential are found phenomenologically. 
They cause the high accuracy of potential models for calculations of the heavy quarko- 
nium masses with no account of spin-dependent splittings, 5m(nL) ~ 30 MeV. Note 
once more, that the scaling properties (|5|), ([?[) and ([11]) are the leading approximations 
in the fitting of the experimentally known masses of the heavy quarkonium levels. 

The accuracy for the predictions of quarkonium wave functions in the framework 
of the potential models is low, for example, it is 5\l/(0)/\I/(0) ~ 30 -f- 50%, since in 
this case the potential behaviour in the border points (r — > and r — > oo) becomes 
essential. One must remember also, that the potential, suitable for the fitting of the 
level masses, can not strictly used for the potential description of the leptonic widths 
(see ref.||). 

In the QCD sum rules, the accuracy of predictions for the heavy quarkonium masses 
is one order of magnitude lower than the accuracy of potential models, 5m sr ~ 
200 -T- 300 MeV. This fact is connected to that the consideration in the QCD sum 
rules takes a finite number of terms in the QCD perturbation theory for the Wilson's 
coefficients and a restricted set of the quark-gluon condensates, so that the results of 
such noncomplete consideration depend on an unphysical parameter, defining a scheme 
of the averaging in the QCD sum rules (the number of moment for a spectral density 
of current correlators or the Borel transform parameter). An additional uncertainty is 
related with a modelling of a nonresonant contribution into the current correlator, i.e. 
with the threshold of hadronic continuum. Such parametric dependences lead to the 
low accuracy of QCD sum rule predictions for the heavy quarkonium masses^. More- 
over, the use of weight functions, defining the averaging scheme and rapidly decreasing 
with the energy rise, causes a suppression of the contribution of higher excitations in 
the quarkonium, so that, as a result, this contribution is neglected, when one tries to 
predict the leptonic constants of the B c meson, for example. 

Recently the QCD sum rule scheme has been offered in papers of refs. || |1J 



so 



this scheme allows one to take into account the spectroscopic characteristics of higher 
S-w&ve excitations. Using the flavour-independence of the S"-wave level density, one 
finds the following regularities. 



1 The QCD sum rule accuracy in calculations of the leptonic constants (f^, fy) is higher (r 
20 -j- 25%), since one uses the heavy quarkonium masses, known experimentally. 
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1) The scaling relation for the leptonic constants / of S-wave levels of the heavy 
quarkonium with the mass M, the reduced quark mass « and the fixed number of the 
quarkonium excitation [ITJ is 

P {M\ 2 . , 

- [-) = const. , (13) 



that is in a good agreement (A/// ~ 5%) with the experimental data [0 on the 



leptonic constants of T, ip and even particles, which are the quarkonia with the 
hidden flavours (QQ), so that 4/x/M ~ 1, and one has [f| 

P 

— = const. , (14) 

independently of the heavy quark flavours in the (QQ) system. Relation (0) essentially 
differs from the scaling law for the leptonic constants of heavy mesons (Qq), containing 
a single heavy quark, where in Heavy Quark Effectice Theory (HQET) [12| one has 

P -M = const. (15) 

Law fllS] ) can be obtained from eq. (|i~3|) in the limit u = m q mQ/(m q + mq) — > m q , 
rriQ ^> m q , M — > mq, so that /i does not depend on the heavy quark flavour. Eq . ([13|) 
gives reasonable estimates for the leptonic constants of B and D mesons 0, if one 



supposes /x ~ 330 MeV [T(J 



2) The scaling relation for the leptonic constants of nS'-levels in the quarkonium is 



(16) 
P 2 ni 



independently of the heavy quark flavours. Eq.(^) is in a good agreement with the 



experimental data on the leptonic constants in the families of ip and T particles [11 
(Af/f < 10%). 

3) Note, that from eq.flTTD one can find the relation for the mass differences of 
nS'-wave levels in the heavy quarkonium 

M n -M 1 Inn 

n > 1 , (17) 



M 2 -M x ln2 



independently of the flavours of heavy quarks in the quarkonium. Eq.(|17D is in a 
good agreement with the experimental data on the masses of ip and T particles |TI 
(5(AM)/AM < 10%), too. 



However, having derived relations ([13|), (jig) , (17), one has used the phenomenolog- 
ical condition, stating the flavour-independence of the heavy quarkonium level density 
and coming from the analysis, made in the framework of the potential models. 
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In the present paper, in the framework of the offered scheme of QCD sum rules, we 
derive eq.(16) and the relation for the S"-wave level density for the heavy quarkonium 



dM n 2A 
-{n = 1) 



dn lnn^ 

where A = m^Qq) — itlq is the difference between the masses of heavy meson and heavy 
quark, n t h is the number of >S-wave levels of the (QQ) heavy quarkonium below the 
threshold of quarkonium decay into the heavy meson pair (QQ) — > (Qq) + (Qq)- In 
the leading order, one has 

A = const. , (19) 



with the accuracy up to power corrections over the inverse mass of heavy quark W2 



(about the role of logarithmic and power corrections, see ref. [0). 

In the leading approximation, stepping from the charmonium to the bottomonium, 
one can neglect a weak logarithmic variation of the number of levels below the thershold, 

\nn th (bb) ~ In n th (cc) . (20) 

From eq.(|TED-(|2"0"D it follows that in the QCD sum rules one can show that 

A 



<T> 



hin^ 



In Section 1 we consider the scheme of QCD sum rules with the account of spec- 
troscopic quantities for the heavy quarkonium and derive relations fllCf ) and (0). In 
Section 2 we make the phenomenological analysis. In the Conclusion the obtained 
results are summarized. 



1 Heavy Quarkonium Sum Rules 

Let us consider the two-point correlator functions of quark currents 

IV(g 2 ) = z|d 4 xe^(0|TJ M (x)jt(0)|0), (21) 
IW) = i J d A xe^(0\TJ 5 (x)4(0)\0) , (22) 



where 



J^x) = Q^Mx) , (23) 
Mx) = Qi(x) l5 Q 2 (x) , (24) 



(25) 



Qi is the spinor field of the heavy quark with i — c,b. 
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Further, write down 

IW) = {-g»„ + q -f-y V {q 2 ) + ^n 5 (g 2 ) , (26) 

where Hy and lis are the vector and scalar correlator functions, respectively. In what 
follows we will consider the vector and pseudoscalar correlators: Hy(q 2 ) and IIp(g 2 ). 
Define the leptonic constants fy and fp 

(Q\J(i(x)\V(\)) = ieWf v M v e ik *, (27) 
(0\J 5 ,{x)\P) = tk^f P e lkx , (28) 

where 

JsM = Qi(x)^ ll Q 2 (x) , (29) 

so that 

(0\Mx)\P) = i hM}p e ik *, (30) 
TOi + m 2 

where \V) and \P) are the state vectors of 1~ and 0~ quarkonia, and A is the vector 
quarkonium polarization, k is 4-momentum of the meson, kp v = Mp v . 

Considering the charmonium (ip, if}' ...) and bottomonium (T, T', T" ...), one can 
easily show that the relation between the width of leptonic decay V — > e + e~ and fy 
has the form 

W - e+e~) = f eW em ^ , (31) 

where is the electric charge of quark i. 

In the region of narrow nonoverlapping resonances, it follows from eqs.(^Tj) - (|30| ) 
that 

-^mU^is) = J2fvnM 2 v J(s-M Vn ), (32) 
^mUp s \s) = EfkMk j- I )2 5(a-MjJ. (33) 

/I n \iit>X i iit>2) 

Thus, for the observed spectral function one has 

-3m<f ( S ) = --smUtp(s) + py P (s, sf) , (34) 

7T 7T 

where p(s, s t h) is the continuum contribution, which is not equal to zero at s > s t h- In 
the following, we will assume that the continuum contribution is equal to the calculated 
perturbative part at s > Sth, so this can lead to an additional parametric dependence 
on s th . 

Moreover, the operator product expansion gives 

U^ CD \q 2 ) = n^ rt V) + C G (<7 2 )(— G 2 ) + GitflimtQiQi) + ... , (35) 

7T 
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where the perturbative contribution Il^ er *^(g 2 ) is labeled, and the nonperturbative 
one is expressed in the form of sum of quark-gluon condensates with the Wilson's 
coefficients, which can be calculated in the QCD perturbative theory. 

In eq.([35]) we have been restricted by the contribution of vacuum expectation values 
for the operators with dimension d — 4. For Co'\q 2 ) one has, for instance, 



CP = TT^T h 'o7.r "' In^- "" 'J - ), (36) 

where 



q 2 /3(3v 2 + 1)(1 - v 2 ) 2 v + 1 _ 9w 4 + Av 2 + 3 
192mim2 q 2 \ 2v 5 v — 1 v 4 



-2 2/ N2 2 1 4 ™l m 2 /o-v 

q = q - (m 1 - m 2 ) , f = 1 ^ — . (37) 

q 2 

The analogous formulae for other Wilson's coefficients can be found in ref.Q. 

In the leading order of QCD perturbation theory it has been found for the imaginary 
part of correlator that [|J 

QmYly^ (s) = a(3ss — s 2 + Qmim 2 s — 2m\s)6(s — (mi + m 2 ) 2 ), (38) 

8lTS 2 

^mUp rt) (s) = ^(s - ( mi - m 2 ) 2 )6(s - ( mi + m 2 ) 2 ) , (39) 

OTIS 

where s = s — mf + m\, s 2 = s 2 — Amis. 

The one-loop contribution into SmII(s) can be included into the consideration (see, 
for example, ref.H). However, we note that the more essential correction is that of 



summing a set over the powers of (ctg/v), where v is defined in eq.(37) and it is a relative 
quark velocity, and as is the QCD interaction constant at the scale of characteristic 
average quark momentum inside the meson ||. In the following consideration, we will 
assume, that the average square of the quark moment is 

(Pq) = MT) , 

where T is the kinetic energy and /x is the reduced mass of the system. Since the 
corresponding virtualities are less than the heavy quark masses, for the numerical 
estimates we use the one- loop expression for the MS scheme "running" coupling in 
QCD with three light flavours nf = 3 

asiifQ)) = (n-2n f /3)\n(p 2 Q )/A 2 QCD ' 

where Aq C d ~ 100 MeV reproduces as{m 2 z ) ~ 0.12 with rif = 5. The Pq value is 
determined in accordance with the result [Q 

(^)=^ 5/3 ((PQ-p'q) 2 ) = ^ 5/3 2(Pq)- 
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In ref. 0] it has been shown that account of the Coulomb-like gluonic interaction 
between the quarks leads to the factor 

m = ^ - — , (40) 

3 v 1- exp( — ^) 

so that the expansion of the F(v) over as/v <C 1 restores, precisely, the one-loop 0(^f) 
correction 

+ ( 41 ) 
3 v 

The additional as correction is related with the hard gluon contribution, it results 
in the factor H, that for the vector state at mi = m 2 = vhq and v — > is equal to 

K Q = 1 - ^ , (42) 



where the scale of the ajj evaluation can be determined in the way, offered in ref.JTJ 
(BLM). Corresponding calculations in the MS scheme lead to the scale e _11//24 m<2 ~ 
0.632mQ |J. The calculation of the if-factor for the pseudoscalar quarkonium with 



m i 7^ m 2 was performed in ref.||15|, where one found 



2«g ( m 2 — mi , vri2 



H P = 1 + ^ [-^- iln^-2 , (43) 

7r \ffi2 + mi mi / 

where one takes at the scale of reduced mass in the system. Result (^) does not 
come to (|42| ) at mi = rri2, and, hence, the hard gluonic corrections to the correlator 
of vector and pseudoscalar currents are different. Since the nonrelativistic QCD of 
heavy quarks, reformulated recently in ref.||16|| and considered in ref.|T7[, results in 



the decoupling of heavy quark spins in the leading order interactions with gluons, the 
mass dependence in the if-factors is determined by the corresponding renormalization 
in full QCD, where it comes from the one- loop calculations of the vertex and quark 



self-energy diagrams. Following ref . , one straightforwardly finds 



H v = l+ 2 ^l( m2 - m 'ln^- 8 -). (44) 
7r \m 2 + mi mi 3/ 

Results (f4~3"D and fl44]) can be obtained from the renormalization factors for vector and 
axial- vector currents of the heavy quark transition Qi — > Q2, considered in ref. ||18||. 
So, one must replace one heavy quark, say Qi, by the antiquark Q\. This leads to the 
subtitutions mi — > —mi, F\ — > H v and Fy — >■ H P with the absolute value prescription 
for the logarith argument. 

The BLM choice of the running coupling scale for mi 7^ m 2 will be considered 
elsewhere. Here we take the scale to be equal to 0.63 v /mim 2 , that gives the estimate 
accuracy, suitable for the current consideration. 
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In accordance with the dispersion relation one has the QCD sum rules, which state 
that, in average, it is true that, at least, at q 2 < 



n J s — q z 



where the necessary subtractions are omitted. QmU^(s) and l^ QCD \q 2 ) are defined 
by eqs.([32t) - (|34] ) and eqs.([35"l) - (|l|), respectively. Eq. fl4~5f) is the base to develop 
the sum rule approach in the forms of the correlator function moments and of the 
Borel transform analysis (see ref.[|J). The truncation of the set in the right hand side 



of eq.p5]) leads to the mentioned unphysical dependence of the fpy values on the 
external parameter of the sum rule scheme. 

Further, let us use the conditions, simplifying the consideration due to the heavy 
quarkonium. 



1.1 Nonperturbative Contribution 

Following refs. 0, [19] , we consider the n-th order g 2 -derivatives of eq.(||) atg 2 = 0. This 
procedure corresponds to the calculation of moments for the spectral densities for the 
current correlators. As was found in refs.0, [19| , the ratio of n-th moment, calculated 
with the account for the gluon condensate, essential for the heavy quarkonia, to the 
n-th moment, calculated in the one-loop approximation of QCD perturbative theory, 
is equal to 

1 + a{n mom )a s - 

~9 (2n mom + 5) (2m Q ) 4 ' ^ ^ 

for the vector states at ni\ = ni2 = vhq. The as term corresponds to the two-loop QCD 
correction. One can see, that the gluon condensate contribution will be essential at 
"large" values of n mom > ni ~ (m^/A) 4 / 3 , where A is the confinement (or condensate) 
scale. For the bottomonium the n\ value is close to 20, so that at n mom < 20 the 
fraction of the gluon condensate contribution is less than 1%, and it rapidly increases 
at n mom > 20 (see refs. 0, |^, [7], |^| , where one can find a manyfold discussion) . Therefore, 
at n mom < rii one can reliably neglect the gluon contributionFj. 



A{n r ' 



1.2 Nonrelativistic Quark Motion 

The nonrelativistic quark motion implies that, in the resonant region, one has, in 
accordance with eq. (|37D , 

v -> . (47) 

2 Unfortunately, the charm quark is not so heavy, and n; ~ 5 so that the region of the used 
approximation is essentially restricted. 
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So, one can easily find that in the leading order 



07T 



so that with account of the Coulomb factor 



F V ) ~ — , 49 
and with the hard gluon correction one obtaines 

Smn'f^^f^ 2 ^. (50) 
For the bottomonium the use of the v — ► limit in expression for the Coulomb 



factor (f49|) in the resonant region is valid with the accuracy less than 5%. As for 
the approximation of the spectral density by expression ([50]), one finds that the ratio 
of moments, calculated with eqs.(|38"D, ( |39"D and ( |50"D for the resonant region s < s t h, 
rapidly tends to 1, so that the maximum deviation from the unit at moment numbers 
n mom ^ x is close to 15%. To decrease the value of the error due to the nonrelativistic 
consideration, M.B.Voloshin || ||, || considers n mom > 8. In the present consideration 
we are satisfied by the accuracy, when the methodics error for the leptonic constants 
is close to 5%, so that we can use nonhard restriction of n mom > 2, when A/// < 8%. 

Note, that in contrast to refs.|5|, ||, ||, we consider the case, when the perturbative 
integrals at s > s t h are compensated by the hadronic continuum contributions in the 
sum rule equations, so that, in practice, we consider the integration over the "resonant" 
region and remember about a parametric dependence on s t h- 

1.3 Phenomenology of Hadronic Contribution 

As for the hadronic part of the correlator, one can write down for the narrow resonance 
contribution 

ds 
s — 

I^V) = f—.Ufpn , M * Pn v, J(*-MgJ, (52) 

J s ~ Q n \ m l + m 2) 

The integrals in eqs.(|51"|)-(|52|) are simply calculated, and this procedure is generally 
used. 

In the presented scheme, let us introduce the function of state number n(s), so that 

n{M 2 k ) = k . (53) 

This definition seems to be reasonable in the resonant region. Then one has, for 
example, that 



rtv s \q 2 ) = Jj^EfvnM'vAs-M^), (51) 



^4 es) ( S ) = sf Vn{s) -y-E^s- Ml k ) . (54) 



A- 
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Further, it is evident that 



l.^0(s-Mi) = ^4-^0(n-k), (55) 



ds ^ ds dn , 



and eg. (|5lD can be rewritten as 



°^V)=/^^^E»(»-*). (56) 



Taking the average value, one finds 

It is evident that, in average, the first derivative of step-like function in the resonant 
region is equal to 

(^e{n-k))~l. (58) 

Thus, in the scheme, that has a rather phenomeno logical extent, suitable for the large 
vrtQ limit, one has 

<nrV)> - /^/vnM^T, (59) 
(4 es) (9 2 )) - I ^ 9/ S2f2pn(s) ^ ^P-. (60) 



s — q 2 (rrii + ^i2) 2 <fe 

To evaluate errors, related with the transformation of the sum over the resonances to 
the integral over the state density, we consider the ratio of n mom -th moments for the 
vector bottomonium states and corresponding continuous approximations for M n and 
f n (see figure [T]). One can conclude that the made transformation gives the stable ratio 
of moments at n mom < 20, and this result is weakly influenced by the variation of the 
continuum threshold s t h'- 5n t h ~ 0.5, and the quark mass threshold is more essential^. 
So, this approximation can lead to the accuracy, close to 10% for the leptonic constants 
estimation. Of course, the absolute value of the integral representation error is related 
with the low numbers of excitations, and its relative contribution rises with the growth 
of the moment number, when the higher excitations become suppressed. 

Note once more, that the performed transformation is purely phenomenological 
representation of the experimental data. 

3 Note, that this procedure can not be considered as a way for the precision evaluation of the heavy 
quark mass, since the corresponding primery accuracy in the continuous description of the level masses 
is low and it can not be better than the value of spin-spin splittings, which are neglected. 
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Figure 1: The ratio of moments for the spectral density of vector resonances 
in the bottomonium, considered as the discret and continuous states, R = 
M dts (n mom )/M con (n mom ). Solid and empry circles correspond to n th = 4 and n t h = 4.5, 
respectively, when m;, = 4.60 ± 0.01 GeV. Dashes show the ratio for m;, pa 4.64 GeV. 
The fitting parameter of the heavy quark kinetic energy equals 0.40 ± 0.03 GeV. 



Further, the g 2 -derivatives of the average resonance contributions in eqs. (p9[) and 
(J) can be rewritten at q 2 = as 



d 



n mom] d q 2n" 

„ mom 



-1)"" 



d n 



n" 



2n m ' 



where 



om ! dq 

b v {n mom ) 
b P (n mom ) 



(4 res) (0)) = b v {n mom ) J 



^(n^(o)) = b P (rC 



ds 



s 2(n mom +l) 

ds 

2(n mom +l) 



( S fvn(s)~^;)\n m ° 

(mi + m 2 ) 



2 (s 2 fk 



dn x 



(61) 
(62) 

(63) 
(64) 



and the averaging is performed with the weight functions, depending on the number of 
the spectral density moment n mom and shown in the right hand sides of eqs. (|6T|) . 
With the n mom growth the by,p quantities will tend to its values at the basic states 



1 

r 



b v,p -»• ~{mi + m 2 )fv uPl 



dM n . ]- J 
-d^ {n = 1) 



(65) 
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where we use M\ »s m\ + 777.2. Of course, the difference between the basic state mass 
and the sum of quark masses becomes essential at the larger values of the moment 
number n mom , where the binding energy of quarks is determined by the value of the 
gluon condensate, say, in addition to the purely perturbative interaction. For the 
bottomonium at n mom ~ 15 — 20, the accuracy of approximation ( |65|) is less than 10%. 

Next, note that in the specified region 2 < n mom < m, one uses the constant value 
of the perturbative density, so that at the largest admissible values of n mom one finds 



n = 1) -T- H PV . (66) 
Mi 7r dn ' \Mj ' v ; 

Moreover, in the considered region of moderate numbers of the spectral density 
moments, the perturbative parts can be rewritten at q 2 = as 

-1 r d n „( nP r^,„, as .„ o „ r ds 



n mom\ ( l q 



2n< 



W n^(0) = g 16/i 2 H Py J g2(n !l +1) . (67) 



Comparing eq.(^) with eqs.(|H]), (0), one can conclude that at 2 < n mom < rii the 
bpy quantities must be independent of n mom with an accuracy less than 15%. As the 
leading approximation, one can stated 

Sm(n ( " r) (s)) = Smn HCD) (s) , (68) 

that gives with account of eqs.(|50|), (|59| ) and ( j60| ) at the physical points s n = M 2 



fpn,Vn = «5 dM n / 4/X 

M n ~ it dn \M n ) 



py 



(69) 



Note, that for the real heavy quarkonia (bb), (be) and (cc), the scale-dependent part in 
the right hand side of eq.(|69"D is approximately constant 



ctsHpy ~ const. , 
and it is practically independent of the total spin of quarks, so that 

hn ^ fpn = fn , (70) 

with the accuracy A/// < 5%. Thus, one can conclude that for the heavy quarko- 
nia the QCD sum rule approximation gives the identity of fp and /y values for the 
pseudoscalar and vector states. 

4 Note, that the calculated imaginary part for the resonance contribution is not assumed to be the 
true physical expression, and it can not be equal to the exact hadronic contribution. As usual in QCD 
sum rules, one assumes only that the calculated expression may be, in some approximation, used for 
the evaluation of the real part of the correlator. In the present paper, one believes that the calculated 
part is close to the average (not exact) hadronic one ( see eq.(pq)). In spite of the absense of explicit 



dependence on the parameters of the QCD sum rule scheme (n mom , s t h), the meaning of eq.(68) is 
strictly defined only in the specified region of moment numbers and after the determination of scheme 
uncertaities, described in the text above. 
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Eq.(m) differs from the ordinary sum rule scheme because it does not explicitly 
contain the parameters, which are external to QCD. The quantity dM n /dn is purely 
phenomenological. It defines the average mass difference between the nearest levels 
with the identical quantum numbers. 

Further, as it has been shown in ref . |2(J , in the region of average distances between 



the heavy quarks in the charmonium and the bottomonium, 

0.1 fm < r < 1 fm , (71) 

the QCD-motivated potentials allow the approximation in the form of logarithmic law 
[0] with the simple scaling properties, so 

(/ " 0(n) , (72) 



n 



i.e. the density of heavy quarkonium states with the given quantum numbers do not 
depend on the heavy quark flavours. 

In ref.0 it has been found, that relation (^) is valid with the accuracy up to 
small logarithmic corrections over the reduced mass of quarkonium, if one makes the 
quantization of 5-wave states for the quarkonium with the Martin potential by the 
Bohr-Sommerfeld procedure. 

Thus, as it has been shown in refs.[|| [K|, for the leptonic constants of S- wave 
quarkonia, the scaling relation takes place 

independently of the heavy quark flavours. 

This conclusion can be drawn with no use of the potential models, since the approx- 
imation, when the difference between the level masses in the heavy quarkonia is flavour- 
independent, is phenomenological observation, leading to the flavour-independence of 
the state density. 

Further, the phenomenology for the quarkonium masses in the framework of the 
potential models leads to the simple scaling relation for the state density ([TT|). 
Note, that from eq.(|ll|) one finds 

M n — Mi = Inn 

M 2 - M 1 In 2 ' 1 ' 

and 

M 2 -M 1 = ^(n = l)ln2. (75) 
dn 

Eq. ([74]) for the differences of n.S-wave level masses of the heavy quarkonium does 
not contain external parameters and it allows direct comparison with the experimental 



data on the masses of particles in the ip and T families [11 
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a(n) 
3 




n 



Figure 2: The experimental values of nS bottomonium (solid dots) and charmonium 
(empty boxes) mass differences a(n) = (M n — Mi)/(M 2 — Mi) and the dependence in 
the present model a(n) = Inn/ In 2. 



Dependence ([74]) and the experimental values for the relations of heavy quarkonium 
masses are presented on figure 0, where one neglects the spin-spin splittings. 

Note, the ^(3770) and ?/>(4040) charmonium states suppose to be the results of the 
3D- and SS'-states mixing, so that the D-wave dominates in the ^(3770) state, and 
the mixing of the 3D- and 3S-wave functions is accompanied by a small shifts of the 
masses, so that we have supposed M 3 = M^( 4040 ) and /f = / 2 (3770) + / 2 (4040). 

As one can see from the figure 0, relation ([F4]) is in a good agreement with the 
experimental data. 

This relation in the hadronic spectra must result in the specific relarions for the 
leptonic constants, if one considers the QCD sum rules. 
Using eqs.([TT|) and (|69D, one gets 



p 

J n\ 
J 112 



Tl2 



(76) 



First, note that eq. ([73]) , relating the leptonic constants of different quarkonia, turns 
out to be certainly valid for the quarkonia with the hidden flavour (cc, bb), where 
4u/M ~ 1 g (see Table [I]). 

Second, eq.(|73|) gives estimates of the leptonic constants for the heavy B and D 
mesons, so these estimates are in a good agreement with the values, obtained in the 
framework of other schemes of the QCD sum rules M. 
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Table 1: The experimental values of leptonic constants (in MeV) for the quarkonia in 
comparison with the estimates of present model. 



quantity 


exp. 


present 


f<t> 


232 ±5 


235 ± 25 


u 


409 ± 13 


408 ± 40 


/t 


714 ± 14 


714 ± 45 



Third, taking a value of the lS'-level leptonic constant as the input one, we have 
calculated the leptonic constants of higher nS'-excitations in the charmonium and the 
bottomonium and found a good agreement with the experimental values (see the thick 
line curvers on figures ^, where we do not take into account the deviation of 4/i/M 
quantity from the unit, however this effect can be essential for the (cc) system. One 
can easily find the modification form of eq. (|76|) 



2 



f 2 Tli M 2 
J n2 1 ni 



(77) 



The leptonic constant values, rescaled under relation (|77|), are presented as the thin 
line curvers on figures [3], |], so that one can see the uncertainty, caused by the small 
excitation approximation. 

These facts show that the offered scheme can be quite reliably applied to the systems 
with the heavy quarks. 

Further, using eq. (|76|) , at q 2 = one can write down 



s th ds „ dn „ r ath , dlnn 



s 



On the other hand, in the leading approximation, one gets 
and, further, 

s th - Si ~ 2MAE , (80) 

where AE = E t h — fnQ — ^q> is the difference between the threshold energies for the 
decay (QQ') — > (Qq) + (Q'q) and the (QQ') pair production. 
In HQET one has 

AE = 2A + 0(l/m Q ) , (81) 

i.e. in the leading approximation one can take AE ~ 2A, being a constant value, 
independent of the heavy quark flavours. 
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Figure 3: The calculated dependence of nS bottomonium leptonic constants and the 
experimental values of /x( n 5)- 



Then one finds 



821 



M ix \M J \nn th 
Comparing eq.(|69l) and eq . (|82[) , one can easily find that in the leading approximation 



(n = l) 



2A 



dn In n t h 

This means, that the average kinetic energy of heavy quarks equals 

A 



a) 



(T) 



lnn 4/l 



Having derived eq.(jS3|), one has assumed, that 

1) the binding energy of quarks in the lS'-state is negligibly small, than the excita- 
tion energy of nS'-levels 

- dM 

E x < A ~ — , (84) 



dn 

2) the excitation energy of levels is small in comparison with the quark masses 



A 



dM 
dn 



< m Q 



J5) 
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Figure 4: The calculated dependence of nS charmonium leptonic constants and the 
experimental values of f^( n s)- 




so that yfs ~ M, 

3) in the leading approximation the hadronic continuum threshold is determined 
by the masses of heavy mesons 

~ M m + M {Qlg) ~ m Q + m Q , + 2A , (86) 

4) the number of states below the threshold is finite and weakly depends on the 
heavy quark flavours. 

Thus, from eq. (|83|) one can conclude that in the framework of the QCD sum rules, 
one gets the universal regularity for the density of S'-wave quarkonium levels, so, in 
the leading approximation, the relation does not depend on the heavy quark flavours 

dM n/ 

— — (n = 1) = const. (87) 



2 Numerical Analysis 

Relation ([83]) is got in the leading approxomation over the inverse mass of heavy quarks, 
when one can neglect the spin-dependent splittings. Therefore, for the A estimate we 
will use the values of S'-level masses of the quarkonia (cc) and (bb). 
One can easily show 

m(n 3 Si) = m(nS) + ^Am(nS) , (88) 

3 

m^So) = m{nS) — -Am(nS) , 
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where Am(nS) is proportional to the leptonic constant squared f% s 0, PD[ , so 
that from the previous Section it follows 

, Am(lS) 

Am n5 = — - . (89) 

n 

From the experimental data one has Am^(lS) = 117 MeV, and, taking into account 
eqs. ([38| [39"1) , one gets 

m^(lS) = 3.068 GeV , 
m v ,(2S) = 3.670 GeV , 

m r (lS) = 9.440 GeV , (90) 
m T (2S) = 10.012 GeV , 

where we have taken into account that 

Am r = Am^p — — - , 
asm 

and as(Y)/cts(ip) ~ 3/4 From eq. (|90|) one has 

(Ma-MO^ ~ 0.602 Gel/, (91) 
(M 2 -M 1 )\ r ^ 0.572 Gel/, 

i.e. in average one has 

(M 2 - Mi) ~ 0.587 ± 0.015 GeV . (92) 

In the T family, where the leading approximation over the inverse heavy quark mass 
must be the most reliable, one has 

n th = 4 . (93) 

Then 

A = (M 2 — Mi) (94) 
Accounting for the variation of nth one gets 

A = 0.59 ±0.07 GeV . (95) 

Note, that the threshold energy of the hadronic continuum can be greater than the 
double mass of heavy meson, since a production of light hadrons (such as 7r mesons) 
can be essential. In the last case the level spacing is rewritten as 

M 2 -Mt = 7 

\n(n th + bn th ) 

2Aln2 SEth 2Aln2 

~ (1 + -T-) + Snth ' 96 

m nth A n th 
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so that the variations of E t h and n t h correlate with each other. The corresponding 
uncertainty of A due to the variation of the hadronic continuum threshold is included 
in the A estimation. 

Estimate (j95|) of the important parameter in HQET is in a good agreement with 



the estimates, made in the QCD sum rules for the heavy mesons [[13 



A ~ 0.57 ±0.07 GeV . (97) 

The accuracy of estimate ( |95|) is within the limits of the accuracy SA ~ 20 MeV, that 
can be achieved, because of the nonperturbative corrections in QCD [ |i~3|j . 
From eq.(^) one finds the estimate 

(/ M (n = 1) ~ 0.85 ± 0.09 GeV , (98) 



dn 



in 



that is slightly greater than the estimates, made in papers of ref . ||22|| , where dM/d 
0.75 GeV was determined in the polinomial interpolation of heavy quarkonium masses. 
Finally, one has 

(T) = 0.43 ± 0.04 GeV . 

Using the 15- level masses of B and D mesons, one can write down the expression 
for the HQET pole mass of heavy quarks 

>»o " lQll (\S) A //; 



2m Q q(lS) 

where /z 2 = 0.4 ±0.1 GeV 2 [23| is the average square of the heavy quark momentum 



inside the heavy meson. Note, that the 0(l/mg) correction to the heavy quark mass 
can be valuable for the charm quark. So, one finds 

m b = 4.69 ± 0.08 GeV , (99) 
m c = 1.30 ±0.12 GeV . (100) 

However, the mass- dependent correction to the HQET pole mass of the charm quark is 
only 1.5cr deviation from the leading approximation 0(1) over the inverse heavy quark 
mass, so that the central value of m c in the equation above can be reasonably taken as 
m c = 1.40 ±0.12, that is more suitable for the considered approximation of a low bind- 
ing in the quarkonium and it agrees with the value, obtained in the consideration of 
the integral representation for the resonance contribution, described above. This effect 
leads also to poor accuracy in the extraction of A value from the charmonium spec- 
troscopy, where the region of n mom , suitable for the made approximations, is strongly 
restricted. 

Next, using the scaling relation, one can get the prediction 

f Bc = 385 ± 25 MeV , (101) 
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that is in agreement with other estimates, performed in QCD sum rules |24| . 
The data on the if) and T particles [1 1]| give 



oq = as 1 



i6«f^ ^\ 2 = 021±a01 

V M 1 J 



3tt 



(102) 



that must be compared with the estimates, corresponding to the described choice of 
scale parameters 



a b = 0.21 ±0.01 
a r = 0.19 ±0.04 



(103) 
(104) 



where the uncertainty in the result for the charm quark is basically related with the 
deviation of (2m c /Mi) 2 factor from the unit as well as with the large contribution of 
0(l/m c ) term in the value of c-quark mass. So, the offered relations are in a good 
agreement with the current data in the limits of the approach accuracy. 



Conclusion 

In the framework of the QCD sum rules and with the use of the phenomenological 
quasiclassical relation for the state density, the expression for the density of S'-wave 
levels of the heavy quarkonium is derived 

dM n 2A 
(n = 1 



dn In n t h 

that in the leading approximation does not depend on the heavy quark flavours. Here, 
A = wi(Qq) — rriQ and n th is the number of nS-levels below the threshold of {QQ') — > 
(Qq) + {Q'q) decay. The analysis of the spectroscopic data on the charmonium and 
bottomonium allows one to do the estimate 



A ~ 0.59 ±0.07 GeV 



that is in a good agreement with the recent estimates from the QCD sum rules for the 
heavy mesons. 

The phenomenology of quarkonium spectra in the framework of the potential models 
leads to the simple scaling expression for the state density, that in the QCD sum rules 
results in the scaling relation for the leptonic constants of excited nS-states 
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